Suppression of the critical temperature for superfluidity near the Mott transition: 

validating a quantum simulator 
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Ultracold atomic gases in optical lattices have proven to be a controllable, tunable and clean 
implementation of strongly interacting quantum many-body systems. An essential prospect for 
such quantum simulators is their ability to map out the phase diagram of fundamental many-body 
model Hamiltonians. However, the results need to be validated first for representative benchmark 
problems via state-of-the-art numerical methods of quantum many-body theory. Here we present the 
first ab-initio comparison between experiments and quantum Monte Carlo simulations for strongly 
interacting Bose gases on a lattice for large systems (up to A'' ~ 3 x 10^ particles) . The comparison 
has enabled us to perform thermometry for the interacting quantum gas and to experimentally 
determine the finite temperature phase diagram for bosonic superfluids in an optical lattice. Our 
results reveal a downshift of the critical temperature as the transition to the Mott insulator is 
approached. 

PACS numbers: 03.75.Hh, 03.75.Lm, 37.10.Jk, 75.40.Mg 



Ultracold bosonic atoms in optical lattices have 
sparked investigations of strongly correlated many-body 
quantum phases with ultracold atoms [I 2 j_ ^that are 
now at the forefront of current research 4 , 5 Jg, 0| ■ For 
increasing interactions between the particles, a bosonic 
superfluid (SF) converts into a Mott insulator (MI), with 
dramatically different properties So far, sev- 

eral of the characteristic ground state properties of the 
systems either in the SF or in the MI have been mea- 
sured 0, [m, [H, [H, Q [S [H, [13, m, m . However, up 
to now, no finite t emp erature phase diagram of the sys- 
tem daiHil,!!!,!!!!!!,!!!! could be determined exper- 
imentally. Furthermore, a series of papers has questioned 
the analysis of the momentum distributions observed in 
the experiments [T2| , arguing that the temperature could 
be higher than anticipated [23, HI] . 

In this paper, we present for the first time a direct 
comparison of an experiment with ultracold bosons in an 
optical lattice with ab-initio finite temperature, quan- 
tum Monte Carlo (QMC) simulations. The simulations 
are performed for realistic trapping potentials and parti- 
cle numbers without free parameters and include impor- 
tant effects such as the finite time-of-flight (TOF) and 
finite imaging resolution [2^. We rely on the sudden 
appearance of narrow interference peaks (or a sudden 
change in the peak width) on top of a broader thermal 
background to detect the onset of long-range phase co- 
herence [SO^I, equivalent to the onset of superfluidity in 
the three-dimensional system. Recently, the question of 
whether a purely normal cloud could also give rise to such 
narrow peaks has been raised theoretically, thus ques- 



tioning the interpretation of this interference pattern as 
a signature for the appearance of a superfluid compo- 
nent d^, HH. So far, this question could not be settled 
due to a lack of an accurate thermometry method in the 
optical lattice potential [s^,!!^- Here, the direct compar- 
ison of the experimental data and the simulations allows 
to determine a temperature for the ultracold lattice gas. 
Theoretically, this temperature can also be computed 
from the measured initial temperature, assuming adia- 
baticity during the loading of the lattice [H, [H, HI, [s^l • 
We find good agreement between the two approaches, 
and conclude that the system is prepared almost adi- 
abatically, up to small non-adiabatic heating effects of 
technical origin which we quantify. For weak interac- 
tions, we observe the sudden appearance of sharp inter- 
ference peaks near a critical temperature Tc, confirming 
the sudden appearance of a condensate as expected for 
weakly interacting gases. For stronger interactions, yet 
below the Mott transition, although the appearance of a 
phase-coherent component with decreasing temperature 
is still clear, a smoother evolution is observed, with broad 
interference peaks also present in the normal fluid (NF) 
phase [1^, [la] . Finally, we establish a method to extract 
the critical temperature Tc for superfluidity in the lattice 
in both regimes. As we approach the Mott transition, a 
downshift of Tc that cannot be accounted for by single- 
particle effects is observed. 

Bose-Hubbard Model with External Confinement. We 
consider a system of ultra-cold bosons in a three- 
dimensional lattice potential with an additional harmonic 
trap. At low enough temperatures, the physics of the sys- 
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tern is restricted to the lowest Bloch band and can be de- 
scribed by a hopping amphtude J and onsite-interaction 
energy U using the Bose- Hubbard model [1, Q 

H = -J^ (alflj +h.c.) 

+ -^^fi^(n^-l)-^fiih^, (1) 

i i 

where a\ (di) creates (annihilates) a particle on the site 
with index i = {ixTiy,iz), fii = a\a.i counts the num- 
ber of particles on site i and denotes the sum 
over next neighbours only. The quantity fj,i = fi — e.^ 
is a difference between the chemical potential fi and 
the confining potential on corresponding site: = 
(m/2) J2a=x,y,z '^a^a^Q' '^i*^ ™ ^^c particle mass, uja = 
2tt X 20— 60 Hz the trap frequency, and da = Aq,/2 the lat- 
tice spacing in the direction a defined by the laser wave- 
length Aq. We thus work in the regime of local density 
approximation, when fii plays the role of the local chem- 
ical potential, the total number of particles N = rii 
being controlled by the global chemical potential /i. In 
the present work, the external confinement is caused by 
a magnetic trap and the Gaussian laser beams creating 
the optical lattice. 

Phase Diagram of the Trapped System The phase di- 
agram for the homogeneous system (i.e. ex,y,z = 0) at 
unity filling is shown in Fig. 1 as a function of interaction 
strength U/J and temperature T. For T — 0, the sys- 
tem undergoes a quantum phase transition (QPT) from 
the SF to the MI phase 8, 9, lOj at the critical interac- 
tion strength (C//J)c = 29.34(2) For interactions 
U/J < {U/J)c, a phase transition between the SF and 
the NF exists at a critical temperature Tc which tends 
to zero as U / J ^ {U/J)c- This behaviour at finite tem- 
peratures can be seen as a generic feature of QPTs, orig- 
inating from a fundamental change in the ground state 
of the system. In this paper, we focus on the transition 
between SF and NF phases and the downshift of Tc. 

In the presence of an external confinement, the system 
becomes inhomogeneous, rendering the notions of global 
'phase diagram' or 'transition' problematic. The decrease 
of the local chemical potential /z^ with distance from the 
trap center results in the coexistence of different phases 
in the trap (see Fig. 2a,b), comphcating the interpreta- 
tion of the experimental data in general. Our strategy to 
avoid this ambiguity is to work with a total number of 
particles N such that the central density stays close to 
unity. In order to achieve this, we determine the chemical 
potential fi = iii corresponding to a central density n = 1 
at T = TciU/J) for a given value of U/J (see Fig. 2c). A 
full QMC study using /ii and including the external trap 
yields the target particle number N which we maintain 
throughout the simulations and experiments. As long as 
U/J < 26.7 (Ta/J > 3), this ensures that the SF starts 
to form first in the center of the trap (see Fig 2a, b) . This 
approach justifies a direct comparison of the critical tem- 
peratures measured for the inhomogeneous system with 
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FIG. 1: Simplified scheme of the finite T phase diagram for a 
single species of bosons in a lattice potential at density n = 1. 
At T = 0, the system undergoes the transition from a SF to 
a MI at the critical interaction strength {U/J)c- For U/J< 
{U/J)c, the SF phase exists up to a critical temperature Tc 
which decreases to zero at the QPC, signaling the drastic 
change in the ground state of the system. The MI phase right 
of the QCP exists strictly speaking only at T = 0. However, 
Mott-like features can be observed at finite temperatures T <^ 
U/kB- The dashed line represents a typical trace in the phase 
diagram along which experimental data and simulations were 
taken in order to determine Tc. 



the ones of the homogeneous system at unity filling which 
are well known from QMC studies [26| . 

Experimental Sequence and QMC Simulations. Our 
experimental sequence starts with a Bose-Einstein con- 
densate (BEG) of ^^Rb atoms produced in a cigar-shaped 
magnetic trap by evaporative cooling. The condensates 
have no discernible thermal fraction and contain a vari- 
able number of TV = 9 x lO"* to 3 x lO^^ atoms. The 
temperature is subsequently varied by exposing the BEG 
to a controlled and calibrated heating sequence, which al- 
lows to access initial temperatures Ti between 20 nK and 
400 nK while keeping N constant (see Methods). Af- 
ter setting the temperature of the gas, we adiabatically 
decompress the magnetic trap towards an expanded, al- 
most spherical trap with radial and axial trap frequen- 
cies of 27r X 18.31(1) Hz and 27r x 11.69(1) Hz, respec- 
tively. Subsequently, the three-dimensional optical lat- 
tice is ramped up within tramp = 300 ms to the final 
depth Vb using an s-shaped ramp [s^. The orthogonal 
retro-reflected laser beams forming the optical standing 
waves have wavelengths of A^; = 765 nm along one di- 
rection and Xy^z = 844 nm along the other two Fi- 
nally, the atoms are released by simultaneously switch- 
ing off the magnetic trap and the lattices and probed 
after ixoF — 15.5 ms using resonant absorption imag- 
ing. This procedure yields the integrated column-density 
n±{x,y) = J nTOF(f) dz which is related to the single- 
particle density matrix in the trap. 

Numerically, the Hamiltonian Eq. ([1]) can be effectively 
simulated by the QMC worm algorithm [4l|, . This is a 
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Vo = 8Er , UIJ=8.^^ , T^°"'= 26.5nK 
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FIG. 2: Phase diagram of the Bose-Hubbard Hamiltonian 
as a function of temperature T/J and chemical potential 
/i/J (a) or density n (b) at the critical interaction strength 
{U/J)c = 29.34(2) (solid line, circles), at U/J = 26 (dotted 
line, squares) and at U/J = 8.5 (dashed line, diamonds). At 
the critical value {^,/ J)c = 11.50(5) (n = 1), the tip of the 
Mott lobe is reached where Tc = 0. In the lower panel (c), we 
plot n/ J corresponding to a density of n = 1 at T = Tc as a 
function of U/J. 



statistically exact method, scaling linearly with the sys- 
tem volume and the inverse temperature. We can deal 
with realistic system sizes (up to ^ 220"^) at the experi- 
mentally relevant temperatures (fcsT < 6 J). For a given 
lattice depth Vq, we calculate the Hamiltonian parame- 
ters J and U from the single particle band structure [i^ . 
The external trap parameters are deduced from the 
lattice depth, the measured laser waists and magnetic 
trap frequencies. Since all parameter values are taken di- 
rectly from the measured experimental control parame- 
ters, this comprises a full ab-initio study. The simulation 
results are translated into integrated column densities af- 
ter TOF, taking into account the finite expansion time as 
well as a finite imaging resolution [2^. The latter is ac- 
complished by convolution of the simulation images with 
the Gaussian point-spread function (5.6 //to root-mean- 
square width) of our imaging system. This value was 
determined experimentally by measuring the autocorre- 
lation function of expanding atom clouds deep in the MI 
regime 

Comparison of QMC simulations and experimental re- 
sults. Two typical sequences of experimental TOF im- 
ages at lattice depths ofVb = SEr iU/J = 8.11) and 
Vo = 11.75 {U/J = 27.5) and for different tempera- 
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FIG. 3: Comparison of experimental and simulated TOF dis- 
tributions: We show the integrated column density n±{x,y) 
represented by the optical density (OD) as obtained from 
the experiment and the QMC simulations for different tem- 
peratures and two lattice depths Vo = 8Er {U/J — 8.11, 
= 2.8x 10^ a-c) and 11.75 £r {U/J = 27.5, N = 0.9x 10^ 
d-f). The simulation results b, e are selected to match the 
corresponding experimental distribution a, d as close as possi- 
ble and the agreement is underlined in the profiles c, f taken 
along the axis denoted in the TOF distributions. The red 
shaded regions in c,f represent the experimental peak-to peak 
fluctuations of the TOF densities, while the red solid line is an 
average over at least three experimental shots. To the exper- 
imental images, we assign the temperatures T/ calculated for 
an adiabatic loading process (see text) where the errors stem 
from the uncertainty in the calibration of the initial tempera- 
ture Ti. The temperatures for the QMC results are exact. 
The comparison at Vo = 8 Er confirms the adiabaticity of the 
loading process, while for Vo = 11.75 we find a small general 
shift in temperature by up to 30% due to heating processes. 



tures are shown in Fig. 3 together with the corresponding 
simulation results. To estimate the temperature of the 
ensemble in the lattice potential, we start from the initial 
temperature Ti measured in the magnetic trap. Using nu- 
merical data for the canonical energies and entropies in 
the initial and final potential, we assign an entropy Si{Ti) 
to a particular experimental run. The final temperature 
is found by inverting Tf = Tf{Si), assuming the initial 
entropy Si{Ti) is conserved during the lattice loading. 

Besides relying on the "adiabatic" temperature T/, an- 
other possibility is to match the experimental profile with 
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that computed from QMC simulations. This yields a 
"matching" temperature T^. Since the temperature used 
for the simulations is exact, this procedure can be seen 
as a direct thermometry for the experiment. In our com- 
parison, the accuracy is limited by the sampling of the 
experimental and simulation data along the temperature 
axis. For Vq = SEr, we find very good agreement be- 
tween the measured and simulated TOF distributions for 
Tf = ry However, for Vq = 11.75 we find the best 
agreement for being slightly higher (up to a maxi- 
mum of 30%) than Tf{Si). We interpret this shift in the 
final temperature as a signal for non-adiabatic heating 
of the ensemble during the loading of the lattice due to 
spontaneous scattering of lattice photons, technical noise 
in the laser setup or the finite ramping times. We can 
estimate the expected heating effect due to the sponta- 
neous scattering of lattice photons and fiuctuations of 
the dipole force (see Methods). For the two final lattice 
depths of 8 Er and 11.75 Er presented in Fig. 3, wc find a 
temperature increase at T = Tc by 3.1(3)% and 30(3)%, 
respectively, which can already quantitatively explain the 
observed shift in the latter case. We note, that at the low- 
est initial temperatures the quantum gas remains well in 
the degenerate regime throughout the parameter range 
investigated. 

For the smallest values of Tf and Tj-, we generally ob- 
serve sharp interference peaks separated by 2hk on top 
of a broader pedestal pattern (Fig. 3). The weight of 
these peaks decreases with temperature until only the 
pedestal is left, which blurs further as the temperature is 
increased even more. In three dimensions, the existence 
of a superfluid requires the existence of a BEC, i.e. long- 
range order [i^. This long-range phase coherence gives 
rise to the narrow peaks in the TOF distributions. We 
therefore use the sudden onset of sharp interference peaks 
(or a sudden change in their width) to identify the onset 
of condensation and superfluidity in the lattice [3Q| • The 
width Wp of the sharp peaks is limited in principle by 
the coherence length Ic of the BEC in the trap according 
to Wp « ht/mlc- In practice, the finite expansion time 
tTOF and imaging resolution do not allow to measure Ic 
accurately when lon g-ra nge order is present. In agree- 
ment with refs. [2^, [3l|, l38|, we find that interference 
peaks can also be observed in the background pattern. 
However, due to the short coherence length in the nor- 
mal phase, these peaks are typically much broader than 
the ones caused by the long-range phase coherent super- 
fluid observed at lower values of U/ J, and vary smoothly 
in temperature throughout the parameter range investi- 
gated here. 

The following evaluation steps have been carried out 
simultaneously for the experimental and the QMC data. 
In general, this evaluation does not require the simulation 
results which we use to verify the method. In order to 
flnd a suitable indicator for the boundary between the 
SF and NF phases, we fit both the experimental and the 



QMC data to an empirical model function 

n±{x, y) = W{x, y) [nc{x, y) + Ubg + nth{x, y)] , (2) 

where the envelope W{x,y) represents the momentum- 
space Wannier function, approximated by a Gaussian. 
Here, we assume that the atomic distribution can be 
split into three contributions, a "condensed component" 
nc{x,y) modeled by Gaussian peaks of width Wp at the 
position of the reciprocal lattice nodes, a constant back- 
ground Tibg modeling an "incoherent component" and a 
"thermal component" nth{x,y) modeling thermal exci- 
tations of the condensate. The latter term is chosen as 
nth{x, y) = (exp[cos(fc^a;) + cos{kyy) - 2\Yf^^ /l'^/'^ . 

In the limit Imax — *■ oo, this corresponds to the analytical 
solution for noninteracting bosons in a lattice potential 
with temperature T = l/ksfi. We truncate the sum 
at Imax = 4 as a compromise to achieve small residuals 
and keeping the numerical effort small at the same time. 
Here, we focus on the fraction of atoms in the Gaussian 
peaks fp and their width Wp alone (46j . 

The fit results for three sets of simulations and experi- 
ments are displayed in Fig. 4 together with the visibility 
of the interference patterns, computed as in ref. [l2|. Al- 
though the fit function Eq.[2]is able to reproduce the TOF 
profiles quite well, there is some ambiguity in the inter- 
pretation of the various components given above. For 
U/J < 20, the thermal and condensed parts are well 
accounted for by the terms nth{x^y) and nc{x,y). Be- 
cause of our finite momentum resolution, the width Wp 
of nc{x, y) is practically constant with temperature up to 
T = T^'^^'^, and the condensed fraction is well evaluated 
from the peak amplitude fp. We find a sudden onset of 
the fraction of atoms in peaks fp as the temperature is 
lowered, where fp increases approximately linearly from 
this point. We fit two lines to the data where one is the 
horizontal axis at /p = 0. The intersection point gives an 
estimate TP*^^*^ of the critical temperature, as was done 
for a harmonically trapped Bose gas in ref. [47*1 . At val- 
ues U / J > 20, the thermal component is not fully cap- 
tured by nth due to interactions between non-condensed 
atoms. To compensate for this, the fit finds non-zero am- 
plitude for ric even in the absence of sharp peaks above 
Tc- However, we find Wp to stay constant until a cer- 
tain temperature T < T^°'^^ above which it suddenly in- 
creases, corresponding to a decreasing coherence length. 
We assign this increased peak width (which can be seen 
directly in Fig. 3) to the thermal component of the lat- 
tice gas. Again, we fit two lines to the data and extract 
a second intersection point at T^^^^^ < TP"'^''. 

We plot TP^^i^ and T^''^*'^ as obtained from the QMC 
simulations and the experimental data as a function of 
U/J in Fig. 5. Here, the experimental data points are 
shifted in temperature to account for the heating during 
the lattice loading (see Methods). Both experimental 
and QMC results consistently exhibit a suppression of 
Tc when the interaction strength is increased towards its 
critical value, as it has been predicted for the homoge- 
neous Bose- Hubbard model [13, H^. Above U/J — 20, as 
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FIG. 4; Fit results for the onset of superfluidity: We plot the 
peak fraction fp (top row), the peak width Wp (center row) 
and the visibility (bottom row) for lattice depths of 8Er (a- 
c, U/J = 8.11), WEr (d-f, U/J = 16.1) and 11.25 & (g-i, 
U/J — 23.7). Both, the results for the experimental data 
(filled circles) and the QMC simulations (open circles) are 
presented. The solid lines in the upper two rows correspond 
to the fits made in order to obtain an estimate for Tc- The 
grey vertical lines indicate Tc for a homogeneous lattice at 
unity filling 



discussed above, this behaviour cannot be extracted from 
the fraction of atoms in sharp peaks due to the inability 
of the fit function to separate the thermal and condensate 
contributions in this parameter range. Instead, it must 
be inferred from the analysis of the peak width, reflect- 
ing the spatial coherence properties of the lattice gas. On 
the other hand, the width alone does not yield enough 
information for lower lattice depth, since the fit routine 
assigns zero amplitude above T^^^^. We thus have to use 
a combination of both methods in order to determine 
Tc, which are found to overlap for intermediate values 
oi U/J (see Fig. 5). The critical temperatures obtained 
from the inhomogeneous system in the experiments and 
simulations are remarkably close to the ones of the homo- 
geneous case with unity filling (solid line in Fig. 5 [1^), 
apart from a small (~ 10 — 15%) systematic shift towards 
lower values. This shift is almost entirely due to applying 
our evaluation procedure to the inhomogeneous density 
cloud in the strongly interacting limit. 

Conclusions and Outlook In conclusion, we present for 
the first time a full quantitative comparison between ex- 
periment and ab-initio quantum Monte Carlo simulations 
for large {N ~ 10^) systems of ultracold bosons in opti- 
cal lattices. Using only experimentally measured param- 
eters as input to the simulations and assuming adiabatic 
loading into the lattice, we find remarkable agreement 
up to U/J ~ 20. Discrepancies in the final tempera- 
ture of up to 30% are observed for deeper lattices, which 
are resolved by accounting for specific heating mecha- 
nisms. The direct comparison of experimental and sim- 
ulated TOF images allows us to perform accurate ther- 
mometry for interacting bosons in an optical lattice. Up 
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FIG. 5: Finite temperature phase diagram and suppression 
of Tc in the lattice. The critical temperature for superflu- 
idity Tc/J as obtained from the measurement (circles) and 
the QMC simulations (diamonds) is plotted versus the inter- 
action strength U/J. Closed symbols mark the values yp<=^'= 
extracted from the peak fraction fp, while the values repre- 
sented by open symbols represent T^"^^^ as obtained from 
the evaluation of the peak width Wp. The experimental data 
has been corrected to account for heating during the lattice 
loading (see text). The solid line is the QMC result for the ho- 
mogeneous Bose-Hubbard model at unity filling, taken from 
ref. [2^. 



to J7/ J ~ 20, we find that for typical parameters the sud- 
den appearance of sharp interference peaks with increas- 
ing temperature yields a reliable measure for the onset of 
superfluidity. For larger interaction strengths, very close 
to the quantum critical region, we observe a smoother 
transition to the NF phase, including broad interference 
features for thermal samples. In this parameter range, 
we find that the change in the width of the interference 
peaks gives more reliable information for the determina- 
tion of the critical temperature Tc- Using these analysis 
techniques, we observe the suppression of Tc upon ap- 
proaching the QCP for the SF-MI transition in both the 
experimental and the simulated data, thus mapping out 
the finite temperature phase diagram of the system. Fur- 
thermore, we find that up to U/J ^ 27.5 the bosonic gas 
remains well in the degenerate regime for the lowest ini- 
tial entropies used, in contrast to the theoretical analysis 
presented in ref. [23 . Our results demonstrate the poten- 
tial of using ultracold atoms in optical lattices to quanti- 
tatively study large-scale condensed-matter physics. The 
direct measurement of the suppression of Tc may further- 
more open the way to approach the region above the QCP 
to experimentally study quantum critical phenomena. 

We would like to thank B. Capogrosso-Sansone, P. N. 
Ma, F. C. Zhang, S. FoUing, H. Moritz, T. EssHnger and 
J. Dalibard for stimulating discussions. This work was 
supported by the DFG, the SNF, the EU (IP SCALA), 
DARPA (OLE program) and AFOSR. The simulations 
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Methods 

A. Controlled heating sequence. 

One important requirement for the presented measure- 
ments is the abihty to change the initial temperature Ti 
of the atomic ensemble without changing the number of 
particles N. We use a one-dimensional optical lattice 
with wavelength Aj, = 844 nm superimposed to the mag- 
netic QUIC trap and perpendicular to its slow axis to 
transfer energy to the ensemble in a controlled way. Af- 
ter evaporative cooling, we slowly ramp this lattice to a 
final value Fheat using a s-shaped ramp [s^, before we 
rapidly pulse it off and on for four times using linear 
ramps of 1 ms length. With the lattice at its high value, 
the created excitations are let to thermalize with the rest 
of the sample over a holdtime of 800 ms, before the lattice 
is ramped down again within 700 ms. After this heating 
sequence, we find no significant reduction of the particle 
number and no residual sloshing could be observed. 

B. Temperature measurement in the magnetic 
trap. 

To measure the temperature Ti of the ensemble af- 
ter the heating sequence, we release the cloud from the 
QUIC-trap and probe its TOF distribution after 18.5 ms 
of free expansion by resonant absorption imaging. The 
integrated column density n±{x, y) is fitted by a bimodal 
distribution (ill 

n{p,z) = ^-I^g,[eMl-p'l<^l-z'l<rl)] 

+nBEc{p,z) (3) 

where p and z are the radial and axial coordinates with 
respect to the orientation of the QUIC-trap in the imag- 
ing plane, gaix) = X^n -^"Z"" the Bose function and 
fiBEciPj ^) is the Thomas- Fermi distribution of the con- 
densed part. From the radial and axial radii Cp and 
az of the thermal component and the trap frequencies 
= 27r X 130.4(2) Hz and = 27r x 16.3(1) Hz, we cal- 
culate the radial and axial temperatures of our sample 
as gil 

Both temperatures differ by maximally a few percent, 
indicating thermalization. To calibrate our heating se- 
quence, we measure Ti as a function of Vheat for each 
total number of particles N used in the experiments. 



C. Adaptation of the tunnel coupling for the 
bichromatic lattice. 

In the experiments we ensure that the redistribution 
of the atoms happens with equal tunnel coupling in all 
three spatial directions of our bichromatic lattice. Based 
on the numerically calculated tunnel couplings, we choose 
the lattice depth in the x-direction to be Vx ~ 1.24 Vy^z + 
0.80 Er throughout the loading of the lattice, where Vy^z 
are the lattice depths of the y- and z-directions. 



D. Estimation of the heating power in the optical 
lattice. 

We estimate the heating rate in the optical lattice po- 
tential following the argument of momentum-diffusion in 
ref. [49] for an atom at rest in a standing wave light 
field. Taking into account the contributions from sponta- 
neous emission and the fiuctuating dipole force, one finds 
a position independent heating rate (E) ~ ErTsmax/'^, 
where F is the lifetime of the excited state of the atom, 
Smax ~ 217^ /A^ is the saturation parameter at an inten- 
sity maximum, A is the detuning of the laser light and f2 
is the coupling Rabi frequency. Given the potential depth 
^sw = ^Asniax/2, this translates into (E) ~ i^rFsc with 
the scattering rate Fgc = Vswr/?iA. 

By comparing QMC simulations to experimental TOF 
images taken after loading of the lattice and a subsequent 
variable holdtime t at a fixed lattice depth Vq, we directly 
measure the increase in the energy per particle versus t. 
From repeated measurements at Vq = 8, 10 an 12 Er, we 
obtain (E) — 0.59(8) Ej-Tsc- When calculating Fgc, we 
take into account the different detunings ^x,y,z of our lat- 
tice beams and the fine splitting of the atomic resonance 
into Di and D2 lines. The reduced measured heating 
power could hint that the excitations created due to the 
heating do not fully thermalize with the rest of the en- 
semble on the timescale of the experiment, possibly due 
to reduced collisions between atoms populating different 
Bloch bands 50]. A quantitative, first principle calcula- 
tion of the heating rate is beyond the scope of this paper, 
where we use the measured value. 

For the measurements presented in the main text, we 
estimate the total heating by assuming the relation (E) = 
0.59 ErTsc to be valid for any lattice depth and making 
use of the almost point symmetric shape of the ramp 
profile [3^. We find a total increase in energy by 9(1) x 
10-^ Er for Vo = 8Er and 12(2) x lO'^ for Vq = 
11.75 Using numerical data for the canonical energy 
in the final lattice potential, we obtain an increase in 
temperature near T = Tc by 3.1(3)% and 30(3)% in the 
two respective cases. 
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